Exactly solvable approximating models for Rabi Hamiltonian dynamics 
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The interaction between an atom and a one mode external driving field is an ubiquitous problem in 
many branches of physics and is often modeled using the Rabi Hamiltonian. In this paper we present 
a series of analytically solvable Hamiltonians that approximate the Rabi Hamiltonian and compare 
our results to the Jaynes-Cummings model which neglects the so-called counter-rotating term in 
the Rabi Hamiltonian. Through a unitary transformation that diagonlizes the Jaynes-Cummings 
model, we transform the counter-rotating term into separate terms representing several different 
physical processes. By keeping only certain terms, we can achieve an excellent approximation to 
the exact dynamics within specified parameter ranges. 



I. INTRODUCTION 

The Rabi Hamiltonian is an elegant model for describ- 
ing the transitions between two electronic states cou- 
pled linearly to a single mode of a harmonic driving field 
within the dipole approximation. Because of its simplic- 
ity in form, it plays an important role in many areas of 
physics from condensed matter physics and biophysics to 
quantum optics P, ■ Given the apparent simplicity of 
this model and its wide range of applicability, it is not 
surprising that various aspects have been studied both 
analytically and numerically |ElllllllIilIllIlllllI| 
Remarkably, exact solutions have not been thus far pre- 
sented except for special cases |j| even though it has 
been suggested that the problem may be solved exactly 
0, 0] The Jaynes-Cummings model is a solvable ap- 
proximation to the spin-boson model that neglects the 
counter-rotating term in the Rabi Hamiltonian |l5l Il6j . 
In general, it provides a reasonable approximation to the 
course-grained dynamics in the limit of weak coupling 
and weak field. For stronger fields and couplings, how- 
ever, the model breaks down. While perturbative treat- 
ments can be used to some extent [ljl [l8| , they give rise 
to fast oscillations and a dependence upon the phase of 
the initial state. Furthermore, it seems rather danger- 
ous to introduce a term as a perturbation which may be 
as strong as terms already present in the unperturbed 
Hamiltonian. 

In this paper, we present an alternative approach for 
including the counter-rotating terms into the Jaynes- 
Cummings model. We do this by transforming the 
counter-rotating term in the Rabi Hamiltonian to the ba- 
sis in which the Jaynes-Cummings term is diagonal and 
then truncating the transformed counter-rotating oper- 
ators to obtain a new series of exactly solvable mod- 
els that are related by various symmetry operations. 
We then compare dynamics of the excited state survival 
probability for our approximating models to the Jaynes- 



Cummings model and to numerically exact solutions of 
the Rabi Hamiltonian and show that our approximate 
models do far better job in capturing both the long time 
decay and fine-structure in both the weak and strong field 
limits. 

The rest of this paper is organized as follows. In SeclTTl 
we obtain and justify the approximating Hamiltonians. 
In Sec. IIIII eigenstates and eigenvalues of these Hamil- 
tonians are found. Excited state survival probability for 
one of the approximating Hamiltonians and its compar- 
ison to the results for the Jaynes-Cummings and Rabi 
Hamiltonains are given in Sec. IIVI 



II. OBTAINING APPROXIMATING 
HAMILTONIANS 

The Rabi Hamiltonian describing interaction of a two 
level atom with a single-mode harmonic field can be writ- 
ten as (h — 1) 



H 



-<j z + va) a + g(o~ + + a )(a + al). 



(1) 



Here cr + and a are spin-flip operators that satisfy 



(J + (J- = 1, 



o o~ = a~a~ = 0, (2) 



a z = 2a + tj~ — 1, and a are the boson creation and 
annihilation operators, and a^a — n is the boson number 
operator. Hamiltonian can be split into two parts as 

H = H JC + V, (3) 
where Hjc is Jaynes-Cummings Hamiltonian 

Hjc — —&z + va)a + g(cr + a + (T~a)) (4) 
and V is the so-called counter-rotating term, 

V = g{(i + a) + <7~a). (5) 

Hjc can be brought to a diagonal form, Hjc, by a suit- 
able unitary transformation, 
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Hjc = UH JC U-\ 



(G) 
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in which U is a unitary operator of the form 19] 

U = e° +A -°~ A \ (7) 

In this paper we restrict our attention to the resonant 
case in which v = uj in Eq.Q. For this case, A has the 
following simple form 



A = 



(8) 



Here we use h for a^a to simplify the notation. The uni- 
tary transformation operator in Eq. Q can be brought 
into the following useful form, 

U = — + K{n){l - a z ) + a + L{h)a - - - a t L(re), (9) 

in which K(h) and L(n) are given by 



K{n) = ( ^L-L ) L(n) = ' 



2^2 



V2n + 2 



(10) 



Here 5(n) is a projection operator on the ground state of 
the field. Unitary transformation generated by U diago- 
nalizes Hjc as follows 

+ | ((<r z + l)Vn+l + (a z - l)Vnj (11) 

in which the eigenstates are given by 

\6i) = \i)\n), \ei) = \])\n) (12) 

where | J.) and | f) are eigenstates of a z with eigenval- 
ues of —1 and +1, while |n) is an eigenstate of n with 
eigenvalue n. 

We now consider how the unitary transformation that 
diagonalizes Hjc transforms the total Hamiltonian 

H = UHU- 1 = UHjcU- 1 + UVU- 1 (13) 

The first term is diagonal and we focus our attention onto 
V = UVU^ 1 . V can be written as a sum of four terms: 



where 



UVU- 1 = V = Vi + % + % + % (14) 

Vt = g{a+F 1 (n)a :i + n-{a)fF 1 {n)) 

V 2 = g(a+a i F 2 {h) + a-F 2 {n)a) 

V 3 = 9(a z (F 3 (n)a 2 + (a^ 2 F 3 (h))) 

V 4 = g(F A {n)a 2 + (a)) 2 F A {h)) (15) 



and Fi , F 2 , F 3 , and F4 are expressed in terms of the K 
and L operators as 

Fi = -L(h)L(h + 2) 



F 2 
F 3 
Fi 



hl + 2V2K(h)) 



(L(n) + L(n + 1)(1 + 2y/2K(h))) 
^= (i(n) - L(n + 1)(1 + 2V2K(n))) 



(16) 



We can distinguish three types of terms in Eq. i|15|) 
based on the physical processes that they describe when 
acting on the \6}^} states in Eq. i|12|) . V\ describes atomic 
excitation or relaxation though absorption or emission of 
three photons. V 2 describes the simultaneous excitation 
of the atom and creation of a photon or simultaneous re- 
laxation of the atom and absorption of a photon. V 3 and 
V4 correspond to creation or annihilation of two photons 
with no net change to the excitation state of the atom. 

The question now becomes whether or not keeping only 
some terms in Eq. I|14[) leads to a solvable model and if 
so, is there a physical justification for keeping only those 
terms? Inspection of Eqs. I|15|) shows that there are two 
obvious cases, 



H\ — Hjc +Vi, H 2 — Hjc + V 2 



(17) 



The reason for their solvability is the same as for the 
Jaynes-Cummings model, viz., there exist pairs of states 
such that the Hamiltonian can induce transitions only 
within each pair. 

To determine if either H\ or H 2 can be used to ap- 
proximate H when describing the system dynamics, we 
will use the same approach that justifies the use of the 
Jaynes-Cummings model as an approximation to the to- 
tal Hamiltonian Thus, we will write H in the inter- 
action picture using Hjc as a free Hamiltonian and then 
analyze oscillatory behavior for different terms. Within 
the interaction picture, H becomes, 



Hj = e l 



iHjaty e —iHjct 



(18) 



Note that any operator that depends only on n and a z 
remains unchanged in the interaction picture. Other op- 
erators that appear in Eq. i|15|l have the following inter- 
action picture form 



( ( T-(at)3) / = (7 -( a t)3 e i« 1 t j 

(trV)/ = a+o)e lW2t . 
(at)? = ( a t)V-<, 



(19) 
(20) 
(21) 



where 



(22) 
(23) 



wi = 2oj — g(y/h + 3 + Vn + 1), 
w 2 = 2uj + g{Vn + 2 + Vfi), 
w 3 = 2uj + g(Vfi - V™ + 2) + |(cr, + 1) 

x (Vn + 3 + Vn + 2-Vn + l- Vfi). (24) 

We can see that oscillation frequencies. u>i are now oper- 
ators. If we expand the states on which these operators 



3 



act in terms of eigenstates of h and a z then we can re- 
place both h and a z with their eigenvalues (n for h and 
±1 for a z ) and Wj's become c- numbers. For states with 
moderate occupation number n we can approximate sums 
of square roots in Eqs. (|22I23() as 

s/n + i + Vn+1 « 2Vn + 1, 

Vn + 2 + v^~2\/n+l. (25) 

Differences of square roots in Eq. (|24|l are of order 1 / 
and the terms involving these differences can be omitted 
if gj ' y/n -C oj. This gives the following approximation for 
the effective frequencies 

wi w 2(w-gVn + l) (26) 
+ gVn + 1) (27) 
w 3 w 2w (28) 

Comparing these effective frequencies, we can see the for 
g > wc h ave 2|(w - gi/n + 1)| < 2uj < 2(w + g\Jn + 1) 
if g-y/n+l < 2w. In this case operators appearing in Eq. 
(|19|l will have the slowest oscillating frequency. Similarly, 
if g < 0, operators in Eq. I|20|) will have slower oscillating 
frequency then operators l|19fl and 121|) if —g\Jn + 1 < 2lu 
is satisfied. Thus, we may expect Hi to give a reasonable 
description of the system dynamics for positive g and H 2 
for negative g for specified ranges of parameters. Even 
though these approximations may be unsatisfactory in 
other regimes, we anticipate that some of the complex 
system dynamics that is present in the Rabi Hamiltonian 
will be manifest in our approximating Hamiltonians. 

We now recall that the sign of g in Eq. (JJJ can always 
be chosen as either negative or positive without the loss 
of generality. This is because there are two unitary trans- 
formation whose action on Hamiltonian JQ| is equivalent 
to changing the sign of g. One is the space inversion 
transformation which changes the sign of a and a 1 and 
leaves a^a invariant. (This transformation is generated 
by exp (iira^a)). Another is the transformation generated 
by exp (if (<r z + 1)) that changes the sign of er + and a~ 
but leaves a z invariant. Thus we can approximate Hamil- 
tonian H by either Hi or H2 depending on our choice of 
sign for g. 

III. EIGENSTATES AND EIGENVALUES OF 
THE APPROXIMATING HAMILTONIANS 

First, we will consider Hamiltonian Hi. Its eigenstates 
and eigenvalues can be found along the same lines as 
for the Jaynes-Cummings model, i.e. by diagonalizing 
suitable two by two matrices. The eigenstates have the 
form 

\4>~) = A n \ t)|n)+S n | |)|n + 3), and 

= B n \ T)|n)-A,| 4)|n + 3). (29) 
Here n > and 

An = J— g- , B n = -7==, (30) 
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FIG. 1: A 2 (n) as a continuous function of n for weak coupling 
(w = = 0.1). 

where 

A n + rj n 
fi n = gVn + 2, 

rj n = 2uj - g(y / n+ 1 + Vn + 3), 

A„ = JlA + Wn- (31) 

Eigenvalues corresponding to eigenstates l|29(l are 
k$ = i((2n+3)w + g(V^+l- Vn + 3)±A„). (32) 

Figffl gives A 2 (n) plotted as a continuous function of n 
in the case weak coupling. It can be seen that for low- 
lying values of n, A 2 in) is close to one indicating that 
in this region eigenstates (|29|l are similar to the Jaynes- 
Cummings eigenstates. 

In addition to eigenstates l|29|l . there are three special 
eigenstates of Hi with eigenvalues 

LO LU ?>UJ 1- 

ko = --, h = --g, k 2 = — - V2g (33) 

All three are also eigenstates of Hjc and are given by 

lxo) = U)|0), |xi) = U)|l), |X2) = |4)|2). 

(34) 

Moving on, we now consider eigenstates and eigenval- 
ues of H 2 . Its eigenstates have the form 

IC) = C n \ |)|n + l) + D„| T)|n + 2), 

\^n) = D„U)|n + l)-C7 n |T)|n + 2) (35) 

Here n > 0. We denote eigenvalues corresponding to 
states \(j>n) by A+. Remarkably, the following relation- 
ship holds between coefficients C„ and A n as well as D n 
and B n viewed as functions of the coupling parameter g, 

C n (g)= An (-g), D n (g)=B n (-g) (36) 
Similarly, we have for eigenvalues 

^(9) = 4(-9)- (37) 



4 



The validity of relations l|36|) and (|37[1 can be easily 
checked if one writes for Hi and H 2 explicit 2x2 ma- 
trices whose diagonalization gives coefhcents in Ea. l|2"§|l 
and Ea. lj3l)l) and eigenvalues for Hi and H 2 . For a given 
n, these matrices only differ by the sign in front of g. 

As in the case of Hamiltonian Hi , there are three ad- 
ditional special eigenstates of H 2 . Two of them have the 
same general form as eigenstates <|35(1 but they do not 
have any simple relation (such as Eqs. (|36I37[0 to the 
special states of Hi. These two eigenstates are 

ir) = cU)|o)+d|t)|i), ie + }^ii>io)- c |T)ii) 

(38) 

Here, c and d are given by 

1 7 

c= — d= , . :, (39) 



where 



7 = 



e = yg 



g + V2(uj + e) ' 



\2gu> 



The corresponding eigenvalues are 

l± = -((J + V2g±2e) 



(40) 



(41) 



The third special state of H 2 is also an eigenstate of Hjc- 
It is given by 



with eigenvalue 



lfo) = |t)|0) 



(42) 



(43) 



With the knowledge of eigenstates and eigenvalues of 
approximating Hamiltonians Hi and H 2 we can calcu- 
late the time evolution of any observable. However, since 
observables of interest and initial states are given in the 
original untransformed picture, it is convenient to remain 
in this picture, in which case the time evolution is deter- 
mined by 



Hi = U~ l HiU, H 2 = U- X H 2 U 



(44) 



Explicit operator forms for Hi and H 2 are given in the 
Appendix. Eigenstates of these operators are obtained 
by acting with U^ 1 on states <|29I34[1 in case of Hi or 
states (|35I38|) in case of H 2 . Eigenstates of Hi are 



16! 



V2 



|> (A n \n + 1) + B n \n + 3)) 



+ -j=\ T} (-£?> + 2) +A»), 



V2 



| l)(B n \n+l) -A n \n + 3)) 



+ — \1)(A n \n + 2)+B n \n)). (45) 



The special states of Hi are given by 

Ixo) = |4)|o), 



1x0 = -^(U)|i)-|T)|o)), 

1X2) = -^(U)|2)-|T>|1». (46) 



For eigenstates of H 2 we have 



1 



\tl>~) = -j=\ i)(C B |n + l)+£> n |n + 3)) 
+-^|T) {D n \n + 2)-C n \n)), 
= -4 1 i)(D n \n + l)-C n \n + 3)) 



V2 



~V2 



\1)(C n \n + 2)+D n \n)). (47) 



The special eigenstates are 

ID - i)(c|0) + -||2)) + -||T)|i), 

m - h>(*)-^|2))--|it>ii>, 
ifo) = iu>ii>+-^it>io>. 



(48) 



We can see that each of the states Q450 and (|47|) is a 
superposition of four eigenstates of operators o~ z and n. 
In contrast, we may recall that eigenstates of Hjq are 
superpositions of only two such states. 



IV. DYNAMICS OF THE ATOMIC SURVIVAL 
PROBABILITY 

We will now consider time evolution of the probability 
P(t) for the atom to be exited if the initial state of the 
system given by | f) |/) = | f, /) where |/) is an arbitrary 
state of the field. P(t) is expressed in terms of (tj z (t)) as 



(49) 



Let us consider the case of Hamiltonian Hi . Straightfor- 
ward calculations using a complete set of eigenstates (J2SJ 
EH give for P(t) 



1 



P(t) = - + 5Re 



n=0 



+2^,1+2 r n 



^A n A n 

* 

+2 Ki ' 

+A n B n+2 F+ +2 *F+e^< + i~<» 



-K+)t 



B n B n+2 F n+2 



-K~)t 



+ -±= [A {0\f)F o -*e^o-^ 



5 



hBoTOFo+V^o-fe)* 
|-Ai<l|/)Ff V^i - *^* 

-B 1 (l|/> J F 1 +*e j(K ^ fe2)t 



Here 

= 

F+ = 



(50) 



\tf) = ^(A n (n\f)-B n (n + 2\f))(51) 

|t,/>- ^(A»(n + 2|/) + B n <n|/»{52) 

(53) 



In order to qualitatively understand time dependence 
of Pit) let us classify contributions from various terms in 
Eq. I|5U|) . All the terms in brackets have the form of time 
dependent exponentials preceded by a factor. Absolute 
values of these factors depend on the initial state of the 
field. The term in the second parentheses is due to the 
overlap of the initial state It?/) with the special states 
(|46|l . Its contribution is negligible for initial states with 
small (0|/), (2|/), and (3|/) components. Oscil- 

lating exponentials that appear in the first parentheses 
can be divided into two groups - those involving differ- 
ences of eigenvalues with the same superscripts and those 
involving differences of eigenvalues with the different su- 
perscripts. 

Let us consider the n dependence of k^ +2 — «+. We will 
again assume that states |/) have not too small average 
n. Using explicit form of eigenvalues given by Eq. Il-Sl'l) 
it can be shown that 



v n+2 



2UJ + 



fJ 



(54) 



Hence, this difference can be approximated by 2u> when 
g/y/n + 1 <§; lj which holds for many couplings and initial 
states of interest. A similar result holds for k~ +2 — k~. 



Thus, the second and third terms in the first parentheses 
in Eq. (| 5 C)|> give primarily oscillating contributions with 
frequency of about 2ui. 



Eigenvalue differences K n+2 — an d K 



n+2 



have 



more complicated n dependences. It can be shown that 
for states with typical n 3> 1 but such that g\Jn + 1 <C oj 



K r, 



2gy/n + 1, 
Auj — 2g\fn 



1. 



(55) 
(56) 



These expressions allow to make connection with the 
standard Jaynes-Cummings model. 

We showed earlier that for weak coupling coefficients 
A n are close to one and, therefore, B n are close to zero for 
moderate values of n (FigQJ. Thus, the first term in the 
first parentheses in Eq. (|50|) is dominant for weak cou- 
pling for states with moderate average n values. Replac- 
ing An — > 1 and B n — > in Eas. (|52l53|l . approximating 
(n + 2\f) with (n\f), using Eq. and neglecting the 

terms in the second parentheses in Eq. (|50|l we obtain the 
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FIG. 2: Excited state survival probability for the Rabi Hamil- 
tonian (RH), the approximating Hamiltonian model (AHM), 
and the Jaynes-Cummings model (JCM) (u = l,g = 0.1) for 
the number initial state of the field with n — 6. 



survival 
model 



robability for the resonant Jaynes-Cummings 



Pit) - ~ 2 



^|<n|/)| 2 cos(2 5v ^TT<) 



(57) 



Using Eq. 1)50(1 , we can compare the survival probabil- 
ity Pahm (t) from our approximating Hamiltonian model 
to Pjcm (t) from the Jaynes-Cummings model as well as 
the exact survival Pjin(t) f° r the Rabi Hamiltonian (ob- 
tained by exact numerical integration of the correspond- 
ing Schrodinger equation). We treat both the Jaynes- 
Cummings model and the approximating Hamiltonian 
model as approximations of the Rabi Hamiltonian. 

Fig. |2 shows survival probabilities P(t) for the three 
models for the weak coupling case and when the initial 
state of the field is the number state with n = 6. All mod- 
els show qualitatively similar behavior of the Rabi type 
oscillations. However, Prh(*) never completely collapses 
(Fig. This effect, although not so pronounced, is vis- 
ible in the approximating Hamiltonian model as well. It 
can also be seen that the Rabi frequency for the Jaynes- 
Cummings model is very slightly larger then the oscilla- 
tion frequency for the Rabi Hamiltonian whereas for the 
approximating Hamiltonian model it is slightly smaller. 

Fig. EH gives survival probabilities for the weak cou- 
pling case with the initial state of the field taken as the 
number state with n = 100. We can see that in this case 
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FIG. 3: Excited state survival probability for the Rabi 
Hamiltonian, the approximating Hamiltonian model, and the 
Jaynes-Cummings model (lu — l,g — 0.1) for the number 
initial state of the field with n = 100. 




. 4 



10 20 30 40 50 60 70 



0)t 



FIG. 4: Excited state survival probability for the Rabi 
Hamiltonian, the approximating Hamiltonian model, and the 
Jaynes-Cummings model (iv = 1, g = 0.1) for the coherent 
initial state of the field with a = 2, ((n) = 4). 



of the strong field both the Jaynes-Cummings model and 
the approximating Hamiltonian model deviate from the 
Rabi Hamiltonian. However, qualitatively, the approxi- 
mating Hamiltonian model gives a better description. As 
the Rabi Hamiltonian, the approximating Hamiltonian 
model shows absence of complete collapse and strong de- 
viation from simple oscillatory behavior. 

Taking the initial state of the field as a coherent leads 
to the results shown on Fig. 0] in the case of the weak 
coupling and weak field ((n) = 4). The approximating 
Hamiltonian model approximates the Rabi Hamiltonian 
better then the Jaynes-Cummings model because it ac- 
count for fast oscillations in P(t) with the frequency of 
about 2ui. However, the intensity of these oscillations is 
weaker compared to the Rabi Hamiltonian. 

In the case of the strong coherent initial field (Fig. 
|5J) The approximating Hamiltonian model again gives a 
better approximation to the Rabi Hamiltonian then the 
Jaynes-Cummings model. Both the Rabi Hamiltonian 
and the approximating Hamiltonian model show almost 
periodic revivals of P(t) with no apparent weakening. In 
contrast to the Jaynes-Cummings model, both models 
do not have a region of nearly constant P(t) before the 
onset of the second group of collapses and revivals that 
is present in the Jaynes-Cummings model. 

It is easy construct explicit expression for P(t) in the 
case of Hamiltonian H 2 using its eigenstates and eigenval- 
ues. One can expect that due to the symmetry properties 



RH 

--- AHM 




FIG. 5: Excited state survival probability for the Rabi 
Hamiltonian, the approximating Hamiltonian model, and the 
Jaynes-Cummings model (o» = 1, g = 0.1) for the coherent 
initial state of the field with a = 8, ((n) = 64). 



given by Eqs. I|36I37|I the following relationship will hold 
between P(t) for Hi and Pit) for H 2 viewed as functions 
of g 



P Hl ig,t) n P H2 (-g,t). 



(58) 



In general, the equality is not exact due contributions 
form the special states for Hamiltonians H± and H% for 
which there are no symmetry relations. We will not pur- 
sue investigation of Ph 2 it) here since for g < it gives 
results that are qualitatively similar to Ph 1 for g > 0. 



Hi 



^=(1 - a z )(l - 5{n))Lin + l)a 2 + h.cA (A.l) 

4V2 / 

Hjc +(- 9 -o+Lin)Lih + 2)a 3 + |<r~ (l + 6{n))a 



il + a z )Lih)a 2 



+ ^=(1 - cr,)(l + <5(n))i(n + l)a 2 + h.c)j (A.2) 



V. CONCLUSIONS 

In this paper we present an approach that allows one 
to add extra terms to the Jaynes-Cummings model of 
an atom in an external field. These additional terms 
add complex oscillatory terms to the survival probability 
which become increasingly important as the field inten- 
sity is increased. By retaining select portions of the full 
counter-rotating term from the original Rabi Hamilto- 
nian we obtain an analytically solvable model that com- 
pares favorably with the numerically exact survival prob- 
abilities from the Rabi Hamiltonian over a wide range of 
parameters even for relatively strong field strengths. 



where operator L{h) is defined in Eq. i|10|) and h.c stands 
for Hermitian conjugate. We can see that in Hamiltoni- 
ans Hi and H 2 , the counter-rotating terms are replaced 
by a number of terms involving various intensity depen- 
dent multi-photon transitions. 
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